REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY) 
24-03-2006 

2.  REPORT  TYPE 

Technical  Paper 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Behavior  of  a  Rocket-Like  Coaxial  Injector  in  an  Acoustic  Field  (PREPRINT) 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

D.W.  Davis  &  B.  Chehroudi  (ERC,  Inc.);  D.G.  Talley  (AFRL/PRSA) 

5d.  PROJECT  NUMBER 

23080533 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/PRSA 

10  E.  Saturn  Blvd. 

Edwards  AFB  CA  93524-7680 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-PR-ED-TP-2  006-089 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S 
ACRONYM(S) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/PRS 

5  Pollux  Drive 

Edwards  AFB  CA  93524-70448 

11.  SPONSOR/MONITOR’S 
NUMBER(S) 

AFRL-PR-ED-TP-2  006-089 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited  (AFRL-ERS-PAS-2006-072) 

13.  SUPPLEMENTARY  NOTES 

Presented  at  the  Institute  for  Liquid  Atomization  and  Spray  Systems  (ILASS)  -  Americas  Meeting,  Toronto,  Canada,  23-26  May  2006. 

14.  ABSTRACT 

A  non-reacting-flow  experimental  investigation  was  undertaken  to  gain  a  better  understanding  of  some  of  the  underlying  physics 
associated  with  the  interaction  of  acoustic  waves  and  a  coaxial-jet  injector  similar  to  those  used  in  cryogenic  liquid  rockets.  Liquid  nitrogen  (the 
round  inner  jet)  and  gaseous  nitrogen  (the  annular  outer  jet)  were  used  under  subcritical,  near  critical,  and  supercritical  chamber  pressures,  with 
and  without  acoustic  excitation.  High-speed  digital  imaging  provided  information  on  the  dynamic  behavior  of  the  jet  under  a  variety  of 
conditions.  It  is  found  that  when  the  jet  is  at  the  pressure  node,  an  externally-imposed  acoustic  field  excites  the  dark-core  of  the  jet  to  a  wavy- 
shaped  structure  consistent  with  the  field’s  characteristics.  Mean  and  root  mean  square  (RMS)  values  of  the  dark-core  length  fluctuations  were 
measured  from  images.  It  is  seen  that  as  the  outer-to-inner-jet  velocity  ratio  increases,  the  RMS  of  the  dark-core  length  fluctuations  decreases 
both  with  and  without  the  existence  of  the  acoustic  field.  It  is  thought  that  a  connection  to  the  rocket  instability  may  be  established  from  these 
data  through  examination  of  the  RMS  values.  It  is  possible  that  decreases  in  the  fluctuation  levels,  shown  to  occur  at  higher  velocity  ratios,  could 
weaken  a  key  feedback  mechanism  for  the  self-excitation  process  that  could  be  driving  combustion  instability  in  rocket  engines.  This  could  offer 
a  possible  explanation  of  the  combustion  stability  improvements  experienced  in  engines  when  a  transition  to  higher  values  of  the  outer-to-inner- 
jet  velocity  ratio  is  made.  Finally,  after  a  careful  review  of  relevant  data  taken  here  and  those  by  others,  there  appears  to  be  a  good  correlation 
between  the  dark-core  length  and  the  momentum  flux  ratio 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE 

OF  ABSTRACT 

OF  PAGES 

PERSON 

Dr.  Douglas  G.  Talley 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

19b.  TELEPHONE  NUMBER 

A 

13 

(include  area  code) 

Unclassified 

Unclassified 

Unclassified 

N/A 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


Unclassified 


Unclassified 


Unclassified 


ILASS  Americas,  19"'  Annual  Conference  on  Liquid  Atomization  and  Spray  Systems,  Toronto,  Canada,  May  2006 


Behavior  of  a  Rocket-Like  Coaxial  Injector  in  an  Acoustic  Field  (Preprint) 


D.  W.  Davis*  and  B.  Chehroudi 
ERC,  Inc.,  Edwards  AFB,  CA 

D.G.  Talley 

Air  Force  Research  Laboratory,  Edwards  AFB,  CA 


Abstract 

A  non-reacting-flow  experimental  investigation  was  undertaken  to  gain  a  better  understanding  of  some  of 
the  underlying  physics  associated  with  the  interaction  of  acoustic  waves  and  a  coaxial-jet  injector  similar  to 
those  used  in  cryogenic  liquid  rockets.  Liquid  nitrogen  (the  round  inner  jet)  and  gaseous  nitrogen  (the  annu¬ 
lar  outer  jet)  were  used  under  subcritical,  near  critical,  and  supercritical  chamber  pressures,  with  and  with¬ 
out  acoustic  excitation.  High-speed  digital  imaging  provided  information  on  the  dynamic  behavior  of  the  jet 
under  a  variety  of  conditions.  It  is  found  that  when  the  jet  is  at  the  pressure  node,  an  externally-imposed 
acoustic  field  excites  the  dark-core  of  the  jet  to  a  wavy-shaped  structure  consistent  with  the  field’s  charac¬ 
teristics.  Mean  and  root  mean  square  (RMS)  values  of  the  dark-core  length  fluctuations  were  measured 
from  images.  It  is  seen  that  as  the  outer-to-inner-jet  velocity  ratio  increases,  the  RMS  of  the  dark-core 
length  fluctuations  decreases  both  with  and  without  the  existence  of  the  acoustic  field.  It  is  thought  that  a 
connection  to  the  rocket  instability  may  be  established  from  these  data  through  examination  of  the  RMS 
values.  It  is  possible  that  decreases  in  the  fluctuation  levels,  shown  to  occur  at  higher  velocity  ratios,  could 
weaken  a  key  feedback  mechanism  for  the  self-excitation  process  that  could  be  driving  combustion  insta¬ 
bility  in  rocket  engines.  This  could  offer  a  possible  explanation  of  the  combustion  stability  improvements 
experienced  in  engines  when  a  transition  to  higher  values  of  the  outer-to-inner-jet  velocity  ratio  is  made. 
Finally,  after  a  careful  review  of  relevant  data  taken  here  and  those  by  others,  there  appears  to  be  a  good 
correlation  between  the  dark-core  length  and  the  momentum  flux  ratio 


Introduction 

Combustion  instabilities  in  liquid  rocket  en¬ 
gines  (LREs)  can  reach  amplitudes  of  100%  of 
chamber  pressure  and  damage  an  engine  in  a 
fraction  of  a  second  [1,2].  The  combustion  insta¬ 
bility  mechanisms  span  from  low-frequency 
mechanisms,  which  is  a  feed-system  coupled 
instability,  to  high-frequency  mechanisms,  which 
depend  on  the  acoustic  modes  of  the  chamber. 
Often  LRE’s  operate  at  chamber  pressures  which 
exceed  the  critical  pressure  of  the  propellants.  In 
such  cases,  while  the  pressure  is  supercritical,  the 
initial  temperature  of  the  propellants  is  often 
subcritical,  but  heats  to  supercritical  tempera¬ 
tures  during  mixing  and  combustion.  The  term 
“transcritical”  has  often  been  used  to  describe 
this  process. 

Many  fundamental  questions  remain  to  be 
answered  regarding  the  physical  behavior  of  the 
propellant  jets  and  other  phenomena  governing 
the  combustion  instability.  One  such  question  is 
how  the  fuel  and  oxidizer  jets  change  behavior 
during  an  unstable  combustion  event.  One  early 
study  by  Heidmann  [3]  investigated  the  behavior 
of  liquid  oxygen  (LOX)  jets  in  a  coaxial-jet-like 
injector  during  a  period  when  combustion  insta¬ 
bility  was  induced.  The  LOX  jet  became  shorter 
once  the  combustor  became  unstable.  Similar 
*  Presently  at  GE  Global  Research, Niskayuna,  NY 
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observations  were  made  by  Miesse  [4]  and  Buf- 
fum  and  Williams  [5]  under  cold  flow  conditions 
when  a  single  round  jet  was  excited  with  an  ex¬ 
ternally-forced  acoustic  field. 

Mixing  of  the  propellant  streams,  if  not  con¬ 
trolling,  is  at  least  intimately  related  to  the  com¬ 
bustion  process.  To  a  first  order  approximation, 
one  measure  of  the  mixing  process  is  the  so- 
called  liquid-core  length,  which  has  been  exten¬ 
sively  studied  by  many  researchers  in  the  past  [6- 
24].  A  compilation  of  experimental  core-length 
correlations,  semi-empirical  theories,  and  other 
data  for  shear  coaxial  injector  studies  involving 
core  length  is  summarized  in  Table  1.  Many  of 
the  correlations  and  semi-empirical  theories  re¬ 
viewed  indicate  scaling  of  the  core  length  with 
outer-jet  to  inner-jet  velocity  ratio  (Vr),  outer-jet 
to  inner-jet  momentum  flux  ratio  ( M ),  density 
ratio  ( p0/pi ),  Reynolds  number  (Re),  and  Weber 
number  (We).  One  difficulty  with  applying  these 
relationships  at  supercritical  pressures  is  that  the 
predicted  core  length  is  very  small  or  zero  in 
magnitude,  because  one  of  the  parameters  in  the 
equation  is  We  and  surface  tension  diminishes 
greatly  or  vanishes  all  together.  To  overcome  the 
difficulty  of  core  length  prediction  associated 
large  We,  one  group  of  researchers  [15-20],  pro¬ 
posed  that  M  can  be  used  to  describe  the  scaling 


of  the  core  length  for  shear-coaxial  injectors,  and 
need  not  include  We  if  the  We  is  sufficiently 
high. 

Chehroudi  and  co-workers  [22-33]  have  in¬ 
vestigated  supercritical  jet  flows  at  higher  Rey¬ 
nolds  numbers  and  turbulent  jets  of  interest  to 
practical  applications  in  propulsion  systems.  The 
single  round-jet  work  by  Chehroudi  et  al.  [26-32] 
was  produced  from  two  different  injectors  made 
from  50  mm  long  sharp-edged  tubes  with  the 
inner  jet  diameters  of  0.254  mm  and  0.508  mm, 
(L/D  of  200  and  100,  respectively).  They  in¬ 
jected  pure  N2  and  02  into  N2,  He,  Ar,  and  mix¬ 
tures  of  CO  and  N2  and  the  arrangement  was 
studied  with  shadowgraphs  and  Raman  imaging. 
It  was  determined  that  the  initial  growth  rate  of  a 
supercritical  jet  was  different  from  that  of  a  sub- 
critical  jet.  Furthermore,  they  quantitatively 
showed  that  the  initial  growth  rate  of  a  super¬ 
critical  jet  behaved  similar  to  variable-density 
gas  jets  [29].  From  this  observation,  a  phenome¬ 
nological  model  of  the  initial  growth  rate  based 
upon  time  scale  arguments  was  proposed  which 
agreed  well  with  the  results  acquired  under  sub- 
critical  and  supercritical  pressures  [29].  It  is  also 
important  to  note  that  Chehroudi  et  al.  [30] 
showed  that  the  growth  rate  measured  from 
shadowgraphs  was  about  twice  as  large  as  the 
growth  rate  of  the  jet  measured  by  Raman  imag¬ 
ing  of  the  jet,  based  on  the  full  width  at  half 
maximum  (FWHM)  jet  thickness  values.  The 
fractal  dimension  of  the  initial  region  of  the  jet 
under  subcritical  and  supercritical  pressures  was 
also  measured  by  Chehroudi  et  al.  [27,  30,  32] 
and  compared  to  the  fractal  dimension  of  other 
liquid  and  gas  jets.  It  was  found  that  the  fractal 
dimension  of  supercritical  jets  was  similar  to  that 
of  gas  jets,  while  the  fractal  dimension  of  the 
subcritical  jets  was  similar  to  that  of  other  liquid 
jets. 

Experimental  Setup 

The  experimental  facility,  shown  in  Fig.  1, 
has  been  described  in  detail  previously  [22-24, 
31,  38],  so  only  a  brief  description  will  be  given 
here.  High-pressure  gaseous  nitrogen  (GN2) 
flows  into  the  system  and  its  flow  rate  is  con¬ 
trolled  with  micrometer  needle  valves.  Tempera¬ 
ture  conditioning  is  accomplished  by  passing  the 
GN2  through  a  shell-and-tube  type  heat  ex¬ 
changer  using  liquid  nitrogen  (LN2)  to  cool  the 
GN2.  The  nitrogen  flows  into  both  the  inner  and 
outer  tubes.  The  temperature  of  the  flow  to  the 
inner  round  jet  of  the  shear-coaxial  injector  is 
considerably  lower  than  that  of  the  outer  annular 
jet,  reaching  the  liquid  (subcritical)  or  densified 


(supercritical)  state.  However,  the  outer  jet  is 
cooled  well  below  the  room  temperature,  but  is 
still  in  gaseous  state.  The  chamber  is  pressurized 
with  GN2  at  ambient  temperature.  For  reference, 
the  critical  pressure  and  temperature  of  N2  is  3.4 
MPa  and  126.2  K,  respectively.  The  pressure 
inside  the  main  chamber  is  maintained  by  adjust¬ 
ing  the  outlet  flow  rate  using  a  triplet  of  valves  to 
provide  an  optimum  level  of  control. 

The  test  article  consists  of  a  chamber  within 
a  chamber.  The  main,  external,  chamber  is  used 
to  create  a  pressurized  environment  and  is  fitted 
with  two  sapphire  windows  for  observation  of 
the  jet  behavior.  The  smaller,  internal  chamber, 
channels  and  focuses  acoustic  waves  onto  the  jet 
at  high  intensity.  The  waves  are  generated  by  a 
high  pressure  acoustic  driver  developed  by 
Hersch  Acoustical  Engineering.  To  obtain  suffi¬ 
ciently  high  amplitudes,  the  driver  must  resonate 
the  normal  modes  of  the  inner  chamber.  There¬ 
fore,  the  frequency  at  which  the  effects  are  stud¬ 
ied  is  limited  to  the  first  two  modes  of  the  inner 
chamber,  being  about  3  kHz  and  5.25  kHz,  re¬ 
spectively.  The  acoustic  wavelength  is  much 
longer  than  the  characteristic  transverse  jet  di¬ 
mensions  at  both  of  these  frequencies  and  the 
conclusions  drawn  were  similar.  Consequently, 
only  the  3  kHz  results  are  reported  here.  The 
injector  is  positioned  at  a  pressure  node  (a  veloc¬ 
ity  anti-node)  to  produce  maximum  velocity 
fluctuations  near  the  jet. 

The  acoustic  driver  was  always  operated  at 
its  maximum  power,  which  was  the  same  for  all 
pressures.  As  a  consequence,  the  acoustic  ampli¬ 
tude  decreased  as  pressure  increased,  due  to  the 
increase  in  density.  The  amplitude  was  about 
183  dB  at  1.49  MPa,  about  180  dB  at  3.5  MPa, 
and  about  178  dB  at  4.9  MPa. 

The  shear  coaxial  injector  used  in  this  work 
(see  Fig.  2)  is  based  on  the  well-characterized 
design  of  the  single-jet  injector  used  in  all  previ¬ 
ous  studies  [24-30]  in  this  apparatus.  The  center- 
post  is  made  from  a  stainless  steel  tube  with  an 
I.D.  (Z)2)  of  0.51  mm  (0.020”)  and  an  O.D.  (D2) 
of  1.59  mm  (0.063”)  with  a  length  of  50.8  mm 
(2.00”).  This  is  used  to  produce  the  inner  jet.  The 
resulting  length  to  inner  diameter  ratio  is  100, 
which  is  sufficient  to  ensure  a  fully-developed 
turbulent  pipe  flow  conditions  at  the  exit.  The 
outer  stainless  steel  tube  creating  the  annular 
passage  had  an  I.D.  (D3)  of  2.42  mm  (0.095”) 
and  an  O.D.  ( D4 )  of  3.18  mm  (0.125”).  This  is 
used  to  produce  the  outer  coaxial  annular  jet.  The 
resulting  mean  gap  width  of  the  annular  passage 
is  0.415  mm  (0.016”),  measured  from  an  image 
taken  at  the  exit  of  the  injector.  The  injector  has 


a  small  bias  of  8%  of  the  mean  gap  width.  The 
recess  of  the  inner  tube  was  0.5  Dh  0.25  mm 
(0.010”)  from  the  end  of  the  outer  tube. 

Because  of  the  proximity  of  the  experimental 
conditions  to  the  critical  point  of  nitrogen,  accu¬ 
rate  measurements  of  temperature  were  neces¬ 
sary  to  obtain  reliable  estimates  of  density  and 
other  quantities  computed  from  the  density.  The 
temperature  profiles  and  the  details  of  the  cali¬ 
bration  were  reported  elsewhere  (Davis  and 
Chehroudi  [22]  and  Davis  et  al.  [23]). 

Shadowgraph  images  of  the  jet  were  taken 
using  either  a  PixelFly  CCD  camera  or  a  Phan¬ 
tom  v5.1  or  v7. 1  CMOS  camera.  The  advantage 
that  the  CMOS  cameras  presented  was  the  ability 
to  produce  high-speed  movies  of  the  jet  for  a 
period  lasting  up  to  several  seconds.  The  framing 
rate  for  the  movies  in  this  work  was  1 8  kHz. 

Results  and  Discussion 

A.  Acoustic  Wave  Interaction  with  Shear- 
Coaxial  Jet.  High-speed  movies  of  the  jet  present 
a  large  amount  of  information  about  the  flow 
being  studied.  The  jets  were  studied  at  three  dif¬ 
ferent  pressures:  a  subcritical  pressure  of  nomi¬ 
nally  1.5  MPa,  a  near-critical  pressure  of  nomi¬ 
nally  3.5  MPa,  and  a  supercritical  pressure  of 
nominally  4.9  MPa.  A  sample  of  ten  consecutive 
images  at  the  subcritical  (rows  one  and  two  from 
top),  near-critical  (rows  three  and  four),  and  su¬ 
percritical  (rows  five  and  six)  chamber  pressures 
are  shown  in  Fig.  3.  Rows  one,  three,  and  five 
are  for  when  the  acoustic  driver  is  off  and  the 
remaining  ones  are  when  it  is  activated  at  ~ 
3kHz.  The  evolution  of  the  jet  in  time  is  from 
left  to  right  in  Fig.  3,  with  the  time  interval  be¬ 
tween  frames  of  55.6  ps.  Prominent  in  all  the 
visualizations  of  the  jet  is  the  existence  of  a  dark 
central  region  which  we  refer  to  as  a  “dark  core”. 
The  dark  core  under  the  unexcited  (i.e.  acoustic 
driver  off)  subcritical  pressures  (Fig.  3,  row  one) 
can  be  approximated  as  a  cylindrical-like  struc¬ 
ture  with  unstable  surface  waves  of  low  ampli¬ 
tude.  However,  upon  increasing  the  pressure  to 
near-critical  and  supercritical  pressures  (rows 
three  and  five)  this  structure  changes  to  a  more 
conical  shape.  The  conical  structure  of  the  dark 
core  has  been  reported  before  for  single-phase 
coaxial  jets  by  Lasheras  and  Hop  finger  [20].  As 
demonstrated  in  images  of  Fig.  3  and  the  other 
visualizations  obtained  in  this  work  (not  shown 
here),  the  conical-shape  structure  was  not  ob¬ 
served  under  two-phase  conditions  (i.e.,  subcriti¬ 
cal  chamber  pressure). 

Excitation  of  the  jet  with  acoustic  driver 
yields  significantly  different  behavior  of  the  dark 


core  compared  to  that  of  the  unexcited  one.  The 
strongest  effect  was  observed  at  subcritical  pres¬ 
sures.  It  is  not  clear  at  this  point  whether  this  was 
due  to  the  subcritical  nature  of  the  flow  or  to  the 
fact  that  the  acoustic  amplitude  was  the  highest 
at  subcritical  pressures. 

It  appears  from  the  movies,  similar  to  those 
used  to  generate  Fig.  3,  that  as  a  portion  of  fluid 
originating  from  the  inner  jet  leaves  the  injector 
tip,  the  momentum  from  the  acoustically-induced 
motion  causes  a  transverse  displacement,  push¬ 
ing  the  core  of  the  jet  into  the  higher  speed  annu¬ 
lar  jet.  The  dense  fluid  from  the  core  is  then  ac¬ 
celerated  by  the  outer  jet.  Upon  reversal  of  the 
acoustic  field,  the  dense  fluid,  initially  from  the 
core  but  now  in  the  high-speed  annular  jet,  ap¬ 
pears  to  maintain  its  transverse  component  of  the 
momentum  imparted  upon  leaving  the  injector 
tip  and  hence  the  dense  fluid  particle  does  not 
reverse  its  direction.  The  dense  fluid  parcel  then 
slows  (both  in  the  axial  and  transverse  direc¬ 
tions)  as  it  arrives  at  the  shear  layer  between  the 
outer  jet  and  the  chamber  fluid  farther  down¬ 
stream,  where  a  “cusp-shaped”  structure  is 
formed  from  the  dense  fluid  originating  in  the 
core  of  the  jet.  Subsequent  mixing  and  heat 
transfer  from  the  outer  jet  to  the  inner  jet  core 
fluid  ultimately  causes  the  fluid  parcel  to  be  in¬ 
distinguishable  from  the  outer  jet  fluid. 

An  investigation  was  also  conducted  to  de¬ 
termine  whether  the  side-to-side  oscillation  of 
the  jet  was  occurring  in  a  3D  helical  mode  or  in  a 
2D  planar  mode.  Visualization  of  the  jet  from  a 
direction  orthogonal  to  the  perspective  of  Fig.  3 
revealed  that  the  jet  oscillation  was  largely  2D  in 
the  direction  of  the  planar  acoustic  waves. 

B.  Dark  Core  Length  Measurement  Tech¬ 
nique.  The  algorithm  for  measuring  the  dark-core 
length  in  this  work  starts  with  an  individual  im¬ 
age  from  which  an  image  histogram  is  computed. 
A  threshold  value  is  then  carefully  defined  to 
identify  the  shape  of  the  dark  core  of  the  jet. 

It  should  be  mentioned  that  confusion  exists 
in  the  literature  especially  when  vague  defini¬ 
tions  are  used  to  define  a  core  length.  The  terms 
potential-core,  potential-cone,  intact-length,  in- 
tact-liquid-length,  and  breakup-length  have  all 
been  used  along  with  various  measurement  tech¬ 
niques.  To  be  clear,  and  to  remove  any  possible 
ambiguity  from  the  data,  the  dark-core  length  is 
defined  here  as  the  connected  dark  fluid  region 
between  the  injector  exit  area  and  the  first  break 
in  the  core  as  defined  by  an  adaptive  threshold¬ 
ing  procedure.  See  Davis  [38]  or  Davis  and  Che¬ 
hroudi  [25]  for  more  details. 


C.  Influence  of  Velocity  Ratio  on  Dark  Core 
Length.  Velocity  ratio  Vr  has  been  a  design  pa¬ 
rameter  for  shear-coaxial  injectors,  particularly, 
as  a  criterion  to  ensure  the  stable  operation  of 
LREs.  For  LOX/FE  engines,  the  design  rule-of- 
thumb  has  been  to  keep  the  velocity  ratio  greater 
than  about  10  to  keep  the  engine  stable,  as  is 
discussed  by  Hulka  and  Hutt  [35].  Although  this 
criterion  has  been  suggested  by  the  experimental 
data,  no  physical  explanation  has  been  provided. 
Related  to  this  is  a  method  to  rate  a  LRE  for 
combustion  instability,  known  as  temperature 
ramping.  Temperature  ramping  is  accomplished 
by  lowering  the  temperature  of  the  H2  while 
maintaining  mass  flow  rate  at  a  constant  value. 
The  lower  the  H2  temperature  is  at  the  onset  of 
the  measured  combustion  instability,  the  greater 
the  stability  margin  that  particular  LRE  is  con¬ 
sidered  to  have.  Temperature  ramping  is  related 
to  the  velocity  ratio  because  one  of  the  effects  of 
lowering  H2  temperature  at  constant  mass  flow 
rate  is  increased  density  and  consequently  lower 
velocity. 

To  determine  the  effects  of  the  outer-jet  tem¬ 
perature  (which  is  GN2  in  this  work)  on  the  coax¬ 
ial  jet,  two  nominal  temperatures  of  ~1 90  K  and 
-140  K  were  studied,  called  “high”  and  “low,” 
respectively.  The  averaged  dark-core  lengths  are 
shown  in  Figs.  4(a)  and  5(a)  as  a  function  of  ve¬ 
locity  ratio  (Vr).  The  RMS  of  the  variations  of 
the  dark-core  length  is  also  shown  in  Figs.  4(b) 
and  5(b).  Figures  4  and  5  present  results  for  con¬ 
ditions  when  the  external  acoustic  field  is  turned 
on  and  off  for  both  nominal  high  and  low  outer- 
jet  temperatures.  According  to  Eroglu  et  al.  [10], 
the  average  of  the  length  measured  from  individ¬ 
ual  images  of  a  set  can  be  regarded  as  the  time 
averaged  value.  It  should  be  noted  that,  for  each 
operating  condition  in  this  work,  length  meas¬ 
urements  were  also  made  from  an  averaged  im¬ 
age  of  a  set.  In  general,  the  length  measured  from 
the  averaged  image  is  slightly  shorter  than  the 
average  of  the  individual  dark-core  values  calcu¬ 
lated  in  a  given  set. 

Evident  in  Figs.  4(a)  and  5(a)  is  the  fact  that 
the  length  of  the  dark  core  decreases  as  the 
chamber  pressure  is  increased.  The  dark  core 
provides  an  indication  of  high-density  regions  of 
the  flow.  At  a  constant  chamber  pressure,  as  Vr  is 
increased,  the  length  of  the  dark  core  decreases 
and  appears  to  approach  a  constant  value.  In  a 
mean  sense,  when  the  dark  core  feels  the  im¬ 
posed  external  acoustic  field,  its  length  is  shorter 
than  or  equal  to  that  when  the  acoustic  driver  is 
turned  off.  At  the  near-critical  and  supercritical 
chamber  pressures,  as  the  Vr  increases,  the  dif¬ 


ference  between  the  lengths  of  the  dark  core, 
measured  with  and  without  the  acoustic  field, 
diminishes.  The  RMS  values  of  the  dark-core 
length  fluctuations,  shown  in  Figs.  4  (b)  and 
5(b),  exhibit  somewhat  similar  trends  to  those 
seen  with  the  mean  values.  It  is  known  that  for  a 
liquid-fueled  rocket,  atomization  and  breakup 
processes,  interactions  between  the  propellant 
jets,  droplet  formation,  and  vaporization  are  all 
affected  by  the  pressure  and,  particularly,  veloc¬ 
ity  fluctuations.  Also,  for  any  chemically¬ 
reacting  system,  the  rate  at  which  energy  is  re¬ 
leased  is  sensitive  to  the  rate  of  change  of  tem¬ 
perature,  density,  pressure,  and,  of  course,  mix¬ 
ture  ratio.  It  is  then  quite  intuitive  to  relate,  in 
some  form,  the  RMS  values  of  the  dark-core 
length  fluctuations  to  mixture  ratio  variations. 
On  the  other  hand,  a  low  RMS  value  can  be  in¬ 
terpreted  as  the  jet’s  inherent  steadiness  (or  in¬ 
sensitivity  to  external  stimuli)  and  vice  versa. 
Examination  of  Figs.  4b  and  5b  clearly  shows 
that  this  property  is  drastically  reduced  as  the 
velocity  ratio  is  increased.  It  is  then  quite  possi¬ 
ble  that  the  observed  improvement  in  combus¬ 
tion  stability  at  higher  values  of  velocity  ratio  is 
a  result  of  the  inability  of  the  jet  to  generate  large 
mass  flow  rate  fluctuations  under  these  condi¬ 
tions,  weakening  a  key  feedback  mechanism  for 
the  self-excitation  process.  In  temperature  ramp¬ 
ing  exercises  for  stability  rating  of  LOX/FE  en¬ 
gines,  the  mass  flow  rate  is  usually  maintained  at 
a  constant  value  [24].  Therefore,  as  the  tempera¬ 
ture  of  the  Fl2  is  decreased  during  a  ramping  epi¬ 
sode,  the  Fl2  becomes  more  dense,  which  de¬ 
creases  the  injector  velocity  ratio.  The  RMS 
plots  shown  here  suggest  that  such  a  decline  in 
this  ratio  amplifies  the  jet’s  inherent  unsteadi¬ 
ness,  providing  a  possible  explanation  for  the 
engine’s  eventual  arrival  into  an  unstable  zone  as 
a  temperature  ramping  test  proceeds. 

D.  Scaling  of  the  Dark  Core  Length.  As 
mentioned  above,  the  core  length  has  in  the  past 
been  scaled  with  many  parameters.  In  some 
works.  We  and  Re  numbers  were  used  to  repre¬ 
sent  velocity  values  in  nondimensionalized 
forms,  since  physical  parameters,  such  as  surface 
tension  and  viscosity,  were  not  varied.  In  other 
works,  based  on  single-phase  results,  velocity 
ratio  (Vr)  or  momentum  flux  ratio  ( M)  were  con¬ 
sidered  to  be  the  scaling  parameters.  The  semi- 
empirical  theory  of  Rehab  et  al.  [  1 7]  produces  an 
equation  suggesting  that  the  length  of  the  core 
scales  with  M'°  5 .  This  equation  was  then  quanti¬ 
tatively  compared  to  single-phase  shear-coaxial 
jet  data,  often  with  equal  densities  [16-18].  Note 
that  for  (inner  and  outer)  jets  of  equal  densities, 


M  ~as  reduces  to  Vr.  The  same  dependence  on  M 
was  also  reported  qualitatively  for  single-phase 
shear-coaxial  jets  of  different  densities  by  Favre- 
Marinet  and  Samano-Shettini  [21].  Additionally, 
Lasheras  et  al.  [19]  stated  the  applicability  of  the 
same  M  ~0J  dependence  for  two-phase  shear- 
coaxial  jets.  However,  they  were  unable  to  make 
core  length  measurements  from  their  images. 

A  plot  of  the  measured  dark-core  length  val¬ 
ues,  of  this  work,  versus  momentum  flux  ratio  is 
shown  in  Fig.  6.  A  clear  distinction  between  the 
subcritical  dark-core  length  (diamond  symbols) 
and  that  for  the  near-critical  and  supercritical 
chamber  pressures  are  seen.  Subcritical  data  in¬ 
dicates  a  much  longer  length  than  at  supercritical 
pressures  for  a  given  momentum  flux  ratio.  It 
should  be  noted  that  the  near-critical  pressure 
data  is  slightly  supercritical,  and  both  the  near- 
critical  and  supercritical  pressure  conditions  pro¬ 
duce  a  single-phase  coaxial  jet.  The  dashed  line 
in  Fig.  6  is  a  least-square  curve  fit  to  the  sub¬ 
critical  data,  and  the  dotted  line  is  a  least  squares 
curve  fit  to  the  near-critical  and  supercritical 
data.  As  indicated  by  the  equations  on  Fig.  6,  the 
single-phase  (i.e.,  near-critical  and  supercritical 
pressures)  data  have  the  same  M  ~n'5  dependence 
form  as  reported  by  others  [16-18,  21].  However, 
the  two-phase  subcritical  data  has  a  weaker  de¬ 
pendence,  M  ,  than  the  single-phase  dark-core 
length.  Other  quantitative  differences  between 
the  subcritical  and  supercritical  cases  have  been 
reported  before.  For  single  round  jets,  Chehroudi 
and  co-workers  [29,  30,  32]  found  that  at  super¬ 
critical  pressures  the  spreading  rate  and  fractal 
dimension  values  were  the  same  as  those  for  a 
gaseous  jet  injected  into  a  gaseous  ambient  with 
different  densities  (i.e.,  variable-density,  single¬ 
phase,  gaseous  jet).  This  is  similar  to  the  present 
dark-core  length  observations.  Under  supercriti¬ 
cal  pressures,  our  coaxial  jet  scales  with  M0'5. 
Therefore,  it  appears  that  this  form  of  the  de¬ 
pendency  on  M  is  not  only  valid  for  gas-gas 
shear-coaxial  jets,  but  for  any  single-phase  shear- 
coaxial  jet. 

Figure  7  was  constructed  with  the  objective 
of  comparing  the  dark-core  length  with  all  of  the 
available  data  that  exists  in  the  literature  for  the 
potential-core  length,  intact-core  length,  and 
breakup-length.  This  figure  represents  all  of  the 
available  data  in  the  literature  concerning  core 
length  spanning  5  orders  of  magnitude  in  mo¬ 
mentum  flux  ratio.  Also,  note  that  as  M  ap¬ 
proaches  zero,  one  reaches  a  limit  defining  a 
single  round  jet  configuration  because  the  outer 
velocity  becomes  zero.  From  data  in  Fig  7,  it 
seems  that  for  M  <  1,  data  points  converge  and 


approach  the  core  length  range  expected  for  sin¬ 
gle  round  jets  reported  by  Chehroudi  et  al.  [36] 
and  Oschwald  et  al.  [34].  The  single-phase  data 
presented  in  Fig.  7  follows  the  dependence  of 
L/Di  =A  AT0'5,  where  the  constant  A  is  between  5 
and  12.  At  M  >  100,  the  experiments  of  Favre- 
Marinet  and  Samano-Shettini  [21]  exhibit  a  re¬ 
circulation  bubble  at  the  end  of  the  core,  and  thus 
the  core  length  decreases.  The  injectors  used  to 
produce  the  two-phase  coaxial  jets  of  Eroglu  et 
al.  [10]  and  Englebert  et  al.  [12]  have  much  lar¬ 
ger  outer  jet  gap  widths  (see  Table  1)  than  what 
is  typical  of  rocket  injectors.  Additionally,  the 
apparatus  of  Eroglu  et  al  [10],  reported  in  Farago 
and  Chigier  [37],  does  not  produce  fully- 
turbulent  inner  jet  until  when  Re  >  104.  The  lack 
of  a  fully-turbulent  inner  jet  and  the  significant 
differences  between  their  injector  and  shear- 
coaxial  ones  used  in  rockets  could  be  the  reason 
why  the  core  length  measured  by  Eroglu  et  al. 

[10]  is  shorter  than  those  observed  in  our  work. 
Englebert  et  al.  [12]  reported  that  the  core  length 
scaled  with  M'0J.  The  core  length  by  Woodward 

[11]  for  the  water  potassium  iodide  solution  with 
helium,  however,  obeys  very  nearly  the  trends 
for  the  subcritical  data  points  (i.e.,  25M°'~  ). 
Considering  that  momentum  flux  ratios  near  10 
are  of  importance  for  LRE,  to  the  best  of  the 
authors’  knowledge,  the  data  for  the  subcritical 
(two-phase)  case  is  the  only  reported  information 
in  the  neighborhood  of  the  M=  10. 

Summary  and  Conclusions 

A  non-reacting  flow  study  of  a  cryogenic 
shear-coaxial  injector  was  conducted  at  pressures 
spanning  subcritical  to  supercritical  values.  The 
flow  from  the  inner  jet  of  the  coaxial  injector 
was  liquid  nitrogen  (or  liquid-like,  if  at  super¬ 
critical  pressures)  and  cold  gaseous  nitrogen 
flowed  from  the  outer  annular  jet,  both  injected 
into  a  chamber  pressurized  with  nitrogen  at  room 
temperature.  The  resulting  coaxial  jet  was  im¬ 
aged  with  a  camera  framing  at  a  rate  of  up  to 
18  kHz.  The  jet  was  excited  transversely  with  a 
high-amplitude  acoustic  driver,  with  the  jet  lo¬ 
cated  at  a  pressure  amplitude  minimum.  The 
following  conclusions  are  offered: 

1 .  The  existence  of  high-amplitude  acous¬ 
tic  waves  alters  the  behavior  of  the  shear-coaxial 
jet.  The  resulting  structure  of  the  jet  exhibits  a 
periodic  shape  corresponding  to  the  transverse- 
velocity  field  created  by  the  acoustic  waves.  The 
periodic  oscillations  imposed  on  the  dark  core  of 
the  jet  is  predominately  in  the  direction  of  the 
transverse  acoustic  velocity.  No  helical  mode  for 
the  jet  was  observed. 


2.  The  root  mean  square  (RMS)  variation 
of  the  dark-core  length  decreases  with  increasing 
velocity  ratio  at  a  given  chamber  pressure  and 
asymptotically  approaches  a  constant  level.  The 
RMS  of  the  dark-core  length  is  greatest  at  sub- 
critical  pressures. 

3.  Previous  research  showed  that  an  epi¬ 
sode  of  so-called  “temperature  ramping”,  used 
for  stability  rating,  could  lead  the  engine  to  an 
unstable  behavior.  Also,  from  other  research,  it  is 
shown  that  coaxial  injectors  with  high  outer-to- 
inner  jet  velocity  ratios  (greater  than  ~10  for 
LOX/H2)  tend  to  be  more  stable.  In  the  current 
work,  the  observation  that  the  RMS  of  the  dark- 
core  length  fluctuations  decreased  at  high  veloc¬ 
ity  ratio  under  both  high  and  low  outer-jet  tem¬ 
peratures  is  considered  to  be  a  potential  explana¬ 
tion  observed  combustion  instability  behavior.  It 
is  possible  that  decreases  in  the  RMS  fluctuation 
levels  could  weaken  a  key  feedback  mechanism 
for  the  self-excitation  process  that  could  be  driv¬ 
ing  combustion  instability  in  rocket  engines. 
This  offers  a  possible  improvement  in  under¬ 
standing  of  the  combustion  instability  in  LRE. 
Ultimately,  this  hypothesis  requires  further  test¬ 
ing  in  a  multi-element,  fired  rocket  experimental 
facility. 

4.  The  quantitative  behavior  of  the  dark- 
core  length  of  the  coaxial  jet  at  near-  and  super¬ 
critical  pressures  follows  a  similar  momentum 
flux  ratio  (M)  dependency  reported  for  the  sin¬ 
gle-phase  shear-coaxial  jets  (i.e.,  12M  "  J).  The 
dark-core  length  for  the  subcritical  chamber 
pressures,  however,  scales  withM-0'2. 

5.  Within  a  range  of  momentum  flux  ratios 
in  between  1  and  10,  the  dark-core  lengths  for 
the  coaxial  jet  under  the  subcritical  (two-phase) 
ambient  pressure  used  here  appear  longer  than 
those  measured  at  the  supercritical  (single-phase) 
condition.  This  range  represents  relevant  values 
for  liquid  rocket  engines,  subcritical  data  define 
a  new  regime  which  was  not  reported  previously. 
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Table  1.  Summary  of  published  operating  conditions,  geometries,  measurement  techniques,  and 
proposed  equations  from  the  literature,  measuring  or  correlating  core  length  for  shear-coaxial  jets. 
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Figure  1.  Picture  of  the  experimental  facility 


Figure  2.  Picture  of  the  shear-coaxial  injector  (left)  and  a  cross-section  drawing  of  the  injector  tip.  The  thermocouple 
used  to  measure  exit-plane  temperature  profiles  is  shown  in  the  drawing. 
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Figure  3.  Consecutive  frames  from  high-speed  shadowgraph  movies  with  the  acoustic  driver  turned  off  (in 
rows  1,  3,  and  5)  and  on  (in  rows  2,  4,  and  6)  at  ~3kHz.  Time  increases  from  left  to  right  with  an  interval  of 
55.6  ms  between  frames.  The  first  two  rows  are  at  a  subcritical  chamber  pressure  (~  1.5  MPa),  the  third  and 
fourth  rows  are  at  a  near-critical  chamber  pressure  (~  3.5  MPa),  and  the  fifth  and  sixth  rows  are  at  a  super¬ 
critical  chamber  pressure  (~  4.9  MPa).  The  acoustic  driver  is  turned  off  for  the  first,  third,  and  fifth  rows 
and  on  for  the  second,  fourth,  and  sixth  at  ~  3  kHz.  The  light  gray  lines  in  the  first  and  second  rows  connect 
fluid  structure  as  they  evolve  in  time. 


Velocity  Ratio 


Figure  4.  Plot  of  the  averaged  dark-core  length  (a) 
and  the  RMS  of  the  length  variations  (b)  normalized 
by  the  inner  diameter.  The  solid  symbols  and  lines 
represent  the  data  for  when  the  acoustic  driver  is  off, 
and  the  hollow  symbols  and  dotted  lines  show  the 
data  when  the  acoustic  driver  was  operated  at  ~ 
3kHz.  The  diamond,  square,  and  up-triangle  symbols 
are  sub-,  near,  and  supercritical  chamber  pressures, 
respectively.  All  cases  are  for  the  high  nominal  outer- 
jet  temperature  of  —190  K.  In  the  inset,  the  words 
sub,  near,  and  super  refer  to  subcritical,  nearcritical, 
and  supercritical  pressure  respectively,  and  the  words 
OFF  and  ON  refer  to  the  acoustic  driver  being  off 
and  on  at  ~3kHz,  respectively. 


Figure  5.  Plot  of  the  averaged  dark-core  length  (a) 
and  the  RMS  of  the  length  variations  (b)  normalized 
by  the  inner  diameter.  The  solid  symbols  and  lines 
represent  data  for  when  the  acoustic  driver  is  off,  and 
the  hollow  symbols  and  dotted  lines  show  results 
when  the  acoustic  driver  is  operated  at  ~  3kHz.  The 
diamond  and  square  symbols  are  for  sub-  and  near- 
critical  chamber  pressures,  respectively.  All  cases  are 
for  the  low  nominal  outer-jet  temperature  of  ~140  K. 
In  the  inset,  the  words  sub  and  near  refer  to  subcriti¬ 
cal  and  nearcritical  pressure  respectively,  and  the 
words  OFF  and  ON  refer  to  the  acoustic  driver  being 
off  and  on  at  ~3kHz,  respectively. 


Momentum  Flux  Ratio 


Figure  6.  Dark-core  length  versus  momentum  flux  ratio.  The  diamond,  circle,  and  up-triangle  symbols  rep¬ 
resent  sub-,  near-,  and  supercritical  chamber  pressure,  respectively.  The  hollow  symbols  are  at  a  high  outer- 
jet  temperature  ( — 190  K)  and  solid  symbols  are  at  a  low  outer-jet  temperature.  The  dashed  line  is  25/M 0  2 
and  the  dotted  line  is  12/M 


100 


6U 

a  _ 

u  - 

-1  =5 

S- 

O  -C 

U  el) 

1 s 

5  -l 

s  0J 

**  u 

S  u 

on  — 
a  « 
u  -a 
-1  a 
<u  S 
o  £ 

u  t 

a  o 
- 

ct 

P 


10 


0.1 


X  * 

X 

4  #  * 

W 

''•n  A 

-J 

25/M"2 

/ 

5  im°-s 

-b 

hS++6- 

■ 

□ 

\  ■■+**. r 

12/M"5 

/ 

■ 

» 

X 

o  Subcritical  P;  High  Outer  T  (*)  o  Nearcritical  P;  High  Outer  T  (*) 

a  Supercritical  P;  High  Outer  T  (*)  ♦  Subcritical  P;  Low  Outer  T  (*) 

•  Nearcritical  P;  Low  Outer  T  (*)  □  Eroglu  et  al.  Re=1456 

■  Eroglu  et  al.  Re=4370  X  Eroglu  et  al.  Re=9328 

x  Favre-Marinet  DR=0.138  Air  ►::<  Favre-Marinet  DR=0.655  Air 

m  Favre-Marinet  DR=0.138  He  m  Favre-Marinet  DR=0.028  He 

Rehab  et  al.  D3/D  1  =  1.37  ■  Rehab  at  al.  D3/D 1  =  2.29 

+  Au  and  Ko  o  Englebert  et  al. 

x  Woodward  KI(aq)  -  N2  •  Woodward  KI(aq)  -  He 

0.01 


0.1 


1  10 

Momentum  Flux  Ratio 


100 


1000 


Figure  7.  Comparison  of  the  present  dark-core  length  measurements  with  all  other  core  length  data  avail¬ 
able  in  the  literature  vs.  momentum  flux  ratio.  Amongst  the  data  reported  by  others,  Eroglu  et  al.,  Englebert 
et  al.  and  Woodward  are  two-phase  flows  and  the  rest  are  single  phase. 


